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A study of various SuSY scenarios is presented in which the Hghtest supersymmetric particle is the gravitino G 
and the next-to-lightest supersymmetric particle is a scalar tau t with lifetimes ranging from seconds to years. 
Gravitinos are interesting dark matter candidates which can be produced in decays of heavier sparticles at the 
International Linear Collider (ILC), but remain undetected in direct searches of astrophysical experiments. We 
investigate the detection and measurement of metastable staus, which may be copiously produced at the ILC 
either directly or via cascade decays. A proper choice of the experimental conditions will allow one to stop large 
samples of f's in the calorimeters of the ILC detector and to study the subsequent decays f — > tG. Detailed 
simulations show that the properties of the stau and the gravitino, such as f mass and lifetime and G mass, 
can be accurately determined at a future ILC and may provide direct access to the gravitational coupling, 
respectively Planck scale. 



1. Introduction 



Supersymmetry (Susy) provides a very attractive scenario to account for the amount of dark matter in the 
universe. If i?-parity is conserved, the hghtest supersymmetric particle (LSP) is stable and thus an ideal dark 
matter candidate. Most frequently discussed is a neutralino LSP. Another interesting candidate is the spin 3/2 
gravitino G, cf. ref. [1], which is part of the supersymmetric particle spectrum in extensions to supergravity, or 
local Susy. The mass of the gravitino is set by the SuSY breaking scale F 

F 

"T-3/2 = ma ~ 



(1) 



where Mp — (SttGn)^^^'^ — 2.4 • 10^^ GcV is the reduced Planck scale and Gn is Newton's constant. The 
masses of the other superparticles depend on the supersymmetry breaking mechanism and on the parameters 
of a specific model. In general TO3/2 is a free parameter in the theory and may extend over a wide range of 
O(eV-TeV). 

If the gravitino is indeed the LSP, it may be produced in cascade decays of heavy SuSY particles. To be 
specific, the next-to-lightest supersymmetric particle (NLSP) is assumed to be a charged slepton, the scalar tau 
f. The dominant decay will then be f ^ tG and since the coupling is gravitational, the lifetime may be very 
long, ranging from seconds to years. The decay-width Tf of the metastable f NLSP is given by 
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This formula will be the basis of the present study. The f decay rate, respectively lifetime tf = P^^, depends 
only on the masses rrif and as well as on the Planck scale Mp ~ no further supersymmetry parameters are 
required. Conversely, the 'supergravity Planck scale' can be expressed in terms of observable quantities 
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A comparison of such a measurement at a particle collider experiment with the value determined in macroscopic 
gravitational experiments, Mp = 2.436(2) • 10^* GeV [2], would be a unique test of supergravity [3]. If the 
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gravitational coupling of the f decay is assumed, the direct measurement of the gravitino mass may be replaced 
by a more accurate determination from the decay- width. Rewriting eq. (2), the f lifetime as a function of the 
stau and gravitino masses can be expressed to a good approximation as 



tf ~ 3.6 • 10^ s 



lOOGeV 



"G 



100 GeV 



(4) 



The cosmological production of gravitino dark matter proceeds essentially via two mechanisms: thermal 
production in 2-body QCD interactions, e.g. gg gG or qg qG, in the reheating phase following inflation [4]; 
and/or late decays of the next-to-lightest sparticle, usually a slcpton or ncutralino, which freezes out as usual 
and decays only after the big bang nucleosynthesis phase, the prediction of which should not be destroyed by 
the hadronic decays [5] . Comparing the calculated gravitino relic density with the observed dark matter density 
puts constraints on the NLSP and gravitino masses and on the parameters for supersymmetric models. Collider 
experiments offer a unique possibility to detect metastable staus and to study the properties of gravitinos, which 
otherwise cannot be observed directly in astrophysical experiments. Some phenomenological investigations to 
detect long-lived staus and to measure their lifetime at the LHC and ILC are discussed in [6-8]. 

In the present study the excellent potential of future ILC experiments and detectors will be explored to observe 
and detect metastable f's, measure precisely their mass and lifetime, and to determine in an independent way 
the gravitino mass. In section 2 the investigated physics scenarios will be presented. Details on experimental 
aspects, detector requirements and event simulation will be given in sections 3 and 4. The experimental analyses 
and case studies are described in section 5. The results are discussed in section 6, followed by conclusions. 



2. Physics scenarios 

Several SuSY benchmark models which have been recently discussed in the literature are chosen for the case 
studies. They represent different supersymmetry breaking mechanisms and a variety of stau and gravitino 
masses leading to different experimental conditions. In general, the models are consistent with cosmological 
constraints, as claimed by the authors. Here they should just serve as valid examples of 'typical' spectra. The 
gravitino masses are somewhat arbitrary. However, most important, the results can be easily transformed to 
other scenarios with different stau and gravitino masses. The masses of the sleptons, light ncutralinos and 
charginos, relevant for the benchmark models under study, are compiled in table I. 

Gauge mediated symmetry breaking (GMSB) usually occurs at rather low scales and a light gravitino is 
naturally the LSP [9]. Typical masses are of order eV to keV which may, however, be extended in the GeV 
range [10]. The GMSB scenario SPS 7 from the Snowmass points [11] is chosen as reference model. It is 
described by the conventional parameters A = 40TeV, Mm = 80TeV, ~ 3, tan/3 = 15 and sign/i = +, 
where A is the universal soft SuSY breaking scale. Mm and Nm arc the messenger scale and index, tan/3 is 
the ratio of the vacuum expectation values of the two Higgs fields, and sign /i is the sign of the higgsino mass 
parameter /i. The gravitino mass is set arbitrarily to wlq = 0.1 GeV. 

In supergravity mediated symmetry breaking (SUGRA) the gravitino mass is a free parameter and of the 
same order as the other sparticle masses [12]. If the NLSP is required to be the charged f, then the common 
scalar mass ttiq in minimal versions has to be small and much lower than the common gaugino mass Mi^2- In the 
scenario FS 600 of ref. [7] the mSUGRA parameters are mp = GeV, M1/2 = 600 GeV, Aq = GeV, tan (3=10 
and sign fi = +, where Aq is the common trilinear coupling. The gravitino mass is chosen to be = 50 GeV. 
The authors of the GDM series [8] require a tighter mSUGRA definition with unified scalar and gravitino masses 
= mQ. The parameters of GDM e are mo = 7713/2 = 20 GeV, A/1/2 = 440 GeV, Aq = 25 GeV, tan/3 = 15 
and sign/x = +. The model GDM C is given by mo = 7713/2 = 100 GeV, M1/2 = 1000 GeV, Aq = 127 GeV, 
tan/3 = 21.5 and sign// = The spectrum of GDM 7/ differs from the previous one essentially only by a lighter 
gravitino, thus mo = m3/2 = 20 GeV, M1/2 = 1000 GeV, Aq = 25 GeV, tan/3 = 23.7 and sign/x = +. 
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Table I: Spectra of sleptons, light neutralinos and charginos [masses in GeV] in SuSY scenarios with a gravitino LSP, 
masses of the first and second generation are identical; the gravitino masses [GeV] used in the case studies are given in 
the last line 
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SPS 7 


FS 600 


GDM e 


GDM C 


GDM 77 


XMSB PI 


xMSB P3 


T\ 


123.4 


219.3 


157.6 


340.2 


322.1 


185.2 


102.5 


T2 


264.9 


406.5 


307.2 


659.2 


652.2 


341.5 


356.9 


V-r 


249.6 


396.4 


290.9 


649.5 


641.5 


327.7 


346.9 


CR 


130.9 


227.2 


175.1 


381.4 


368.5 


192.9 


109.5 


Cl 


262.8 


405.6 


303.0 


662.7 


655.3 


340.1 


357.4 
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247.4 


X? 
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243.0 
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426.5 


203.6 


189.2 


xl 


277.9 


469.6 


338.2 


801.9 


801.5 


385.5 


263.8 


xf 


275.5 


469.9 


338.0 


801.9 


801.4 


388.2 


251.5 


G 


0.1 


50 


20 


100 


20 


50 


25 



In gaugino mediated symmetry breaking (xMSB) the gravitino may be the LSP with a mass comparable to the 
other sparticles [13]. In the scenarios proposed by [14] a gravitino LSP can be accompanied by a stau, sneutrino 
or ncutralino NLSP, where the latter options are essentially unobservablc. The spectra of the f NLSP models 
arc determined by the parameters M1/2 = 500 GeV, tan/3 ~ 10, sign/i = + and the soft Higgs masses, chosen 
as ml = ml = OTeV^ for point yMSB PI, and m? = OTeV^ to? = 0.5 TeV^ for point vMSB P3. An 
interesting feature of this particular gaugino mediation scenario is that a naive dimensional analysis provides a 
lower bound on the gravitino mass of m3/2 10 GeV [15]. In the present case studies values of m^ — 50 GeV 
and TOq = 25 GeV are chosen for PI and P3, respectively. 

The sparticle spectra of table I are calculated from the original set of parameters using the programs 
SPheno [16] for models 1, 2 and SuSpect [17] for models 3-7. 

3. Stau detection & measurement principles 

The programme to determine the properties of the metastable stau and to measure, independently, the mass 
of the gravitino consists of several parts: 

- identify f by the characteristic heavy ionisation —dE/dx cx l//3^ in the TPC; 

- determine f mass from two-body kinematics of e+e^ — > fifi production; 

- follow low momentum f candidate until it stops inside detector, record f location and time stamp to; 

- trigger decay f tG at ttrig uncorrelated to beam collision; 

- associate vertex to previously recorded stopping point, get f lifetime tf — ttrig — to] 

- measure r recoil spectra in calorimeter, get gravitino mass TOq. 

The requirements on the detector performance and operation will be discussed in the following. 

The detector is taken from the Tesla tdr [18]; a cross section through one quadrant is displayed in figure 1. 
This concept serves as a baseline for a future LDC detector [19]. In particular the calorimeter layout and 
depth may be taken as representative for various detector designs discussed in the ILC community. The main 
characteristics of the detector, relevant to the present study, are: 

- A TPC with excellent tracking performance, i.e. momentum resolution S{l/pt) < 2 • lO^^GeV^^ and 
dE/dx resolution < 5 %. 
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Figure 1: View of one quadrant of the Tesla detector [18], length units in mm; the amount of material integrated 
perpendicular and along the beam direction is indicated by slightly larger (red) symbols as R [gcm^'^] 



- A highly segmented hadronic calorimeter (HCAL). 40 samplings with 3x3 cm^ lateral readout cells, leading 
to energy resolutions of 6Ehad = 0.5^/E/GcV © 0.04 S for hadrons and SE^m = Q.2^/eJGcV © 0.02 E 
for electromagnetic showers. 

- An instrumented iron flux return yoke equipped with several layers for muon detection as well as for 
calorimetric measurements of hadron showers. 

The amount of absorber material integrated along the radial and longitudinal directions, R [g/cm^], can be 
read off the detector figure 1. The range of a heavy, non-strongly interacting charged particle of mass m, is 
conveniently expressed by R/m [gem" -^GeV"^] and can be parametrised as a function of the scaled momentum 

p/m ~ 

^ogiaiR/m) = ci + C2logio(p/m) , (5) 

with ci = 2.087 and C2 = 3.227 for steel [20]. Values of (3j below which a heavy f can be stopped in the hadron 
calorimeter or instrumented iron yoke depend on the location and angle of incidence. Typical ranges are listed 
in table II. 



4. Event generation 



Events are generated with the program Pythia 6.3 [21] which includes initial and final state QED radiation 
as well as beamstrahlung [22]. The detector simulation is based on the detector proposed in the Tesla tdr [18] 
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Table II: Acceptance ranges of maximal /37 at which heavy r's can be absorbed in the HCAL or instrumented iron yoke 
of the ILC detector 

ruf p-y (HCAL) p-y (Yoke) 

125 GeV 0.41 - 0.46 0.52 - 0.59 
250 GeV 0.33 - 0.37 0.42 - 0.48 
375 GeV 0.29 - 0.33 0.37- 0.41 



and implemented in the Monte Carlo program Simdet 4.02 [23]. 

The analysis of SuSY scenarios with a gravitino LSP is quite different from conventional supersymmetric 
models. The metastable f NLSP will be detected before it leaves or will be absorbed in the detector. Thus, 
there is a very clean signature without missing particles (except neutrinos from decays): the observed particle 
momenta are balanced but their moduli don't sum up to the cms energy 



and Pi < \fs . 



These features are very distinct from Standard Model background which can be efficiently rejected. They also 
allow the sparticle production and decay chains to be reconstructed from the complete event kinematics. Each 
Susy event contains two f 's which can be readily identified as highly ionising tracks and their passage through 
the detector can be accurately followed. The location of stopping f's may be determined within a volume of a 
few cm'^. 

The decay f tG is triggered by an isolated, high energy hadronic or electromagnetic cluster in the HCAL 
(threshold Eh, em > 10 GeV), by a hadronic shower in the yoke {Eh > 10 GcV) or by an energetic /i originating in 
the HCAL or yoke {Efj_ > 10 GeV), occurring at any time not correlated to beam collisions. The main background 
is expected to come from cosmic rays and can be rejected by excluding decay vertices in the outermost detector 
layers or signals being initiated by muons from outside. Further excellent background discrimination is provided 
by requiring the decay vertex to coincide with a previously recorded r stopping point. The t decay modes 
contributing to the analysis are the leptonic 3-body decays r M^m^^t (17.4%) and r eVeV^ (17.8%), and 
the hadronic decays r — * ni^r (11.1%), t pvr n^n^i^r (25.4%) and t iirvr tt^tt^tt^ Vr + 7r*7r°7r°i/,- 
(19.4%). 

The T recoil energy gives access to the gravitino mass 

i^r = ^ ^1 - for < m? . (6) 

The sensitivity to low gravitino masses decreases rapidly, e.g. for ~ 0.1 and rrif precisely known the 
energy E^., which is also the maximum energy of the observable r decay products, has to be measured with a 
precision well below one percent. 



5. Experimental analyses - case studies 

A big advantage of the ILC is that the centre of mass energy can be adjusted in order to control the sparticle 
production and the number of staus trapped in the detector. A range-out in the calorimeter is preferred for 
measuring the kinematics of f decays and thus the gravitino mass, while all staus stopped somewhere in the 
detector can be used for determining the lifetime. Low momentum f 's with a suitable /?7 factor can be produced 
either directly or via cascade decays from light sleptons, e.g. e^j erfi, or neutralinos, Xi ~* 't^i- However, 
all these processes — rifi, eueji, P^rP^r and XiXi — rise only slowly above kinematic threshold with cross 
sections a <x. 0^, thus providing relatively low rates. Obviously, choosing fifi pair production as the only source 
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is not optimal. A more efficient reaction, if kinematically accessible, is the associated selectron production 
e^e~ — *■ eflCL erfi ex^j which increases much faster as cr oc /3 near threshold [24]. Note, such a threshold 
behaviour also applies to e~e^ eRCR production [25]. 

For the present studies the collider energies are not always optimised. Instead, the often proposed 'standard' 
energies of 500 GeV (models 2, 3, 7) and 800 GeV (models 4, 5) are assumed, as well as some special settings 
around 420 GeV (models 1, 6). The integrated luminosities arc moderate and chosen such as to enable a 
reasonable measurement of the f decay spectra. The results obtained may be easily scaled to apply to different 
running conditions, event statistics and detector efficiencies. 

The experimental analysis and techniques will be discussed representatively in some detail for the GDM e 
scenario. For the other benchmark models the main features and results will be presented in a more compact 
form. All results are summarised in table III. 

5.1. mSUGRA scenario GDM e 

The assumptions for the case study of the GDM e scenario consist of the obscrvables m.f ~ 157.6 GeV, 
tf = 2.6 • 10^ s, ttiq = 20 GeV and the experimental conditions of ^/s = 500 GeV centre of mass energy, 
C = 100 fb^^ integrated luminosity, and asuSY = 300 fb as inclusive f production cross section. 

The scaled momentum distribution of f's produced in various reactions is shown in figure 2 a. The 
majority of particles, coming from diagonal slepton and neutralino pairs and peaking around /?7 ~ 1, will leave 
the detector. One observes, however, a second peak at low ^ 0.5 which stem from cascade decays of ejiij^ 
production and which may be stopped inside the detector, see ranges in table II. One expects iV^^^^^ = 4100 
and N^°^'^ = 1850 f's trapped in the calorimeter and yoke, respectively. 

The stau mass measurement is based on the kinematics of prolific pair production e+e^ — > fifi, see magenta 
curve in figure 2 a, to be identified as a pair of coUinear, non- interacting particles with momenta < ~ Ef, 

i.e. well below the beam energy. Determining the mean value of the f momentum with an accuracy of 
(pf) = 192.4 ± 0.2 GeV leads to a precise f mass of 

mf = 157.6 ±0.2 GeV . 

The stau lifetime measurement is based on the decays of f's which have been stopped in the detector. Requir- 
ing an isolated energetic cluster above a threshold of 10 GeV or a muon above 10 GeV originating somewhere 
inside the sensitive fiducial volume of the calorimeter or yoke, results in the decay time distribution shown in 
figure 2 c. A fit to the spectrum gives a f lifetime of 

tf = (2.6 ±0.05) • lOS , 

corresponding to roughly one month. The relation between stau lifetime and gravitino mass as a function of 
stau mass is graphically presented in figure 2 d. For the present condition the lifetime increases as t^ ^ 
and rises significantly faster for ttIq ^ 0.25 m^, reaching 10^ s for a 75 GeV gravitino. Comparing with other 
physics scenarios one also observes that the lifetime can get very long for light f's and small mass differences 
and, for a fixed gravitino mass, becomes shorter the heavier the stau. 

A direct gravitino mass measurement can be performed by exploiting the r recoil of the decay f — > tG, see 
eq. (6) and discussion in sect. 4. The upper endpoints of the energy spectra which coincide with the primary t 
energy Et — 77.5 GeV, are directly related to the masses involved. The leptonic 3-body decay r evv is not 
very useful due to the soft spectrum peaking at low values. The hadronic 2-body decay t ~* ttv produces a flat 
spectrum, where the upper edge is, however, strongly diluted due to resolution effects and limited statistics. 
Well defined upper edges are provided by the hadronic decays to heavier final states t pv and r — > mnriy. 
The energy distribution of both decay modes, defined as 'r jets', is shown in figure 2 b. In order to illustrate the 
sensitivity to the gravitino mass simulations assuming the nominal value of — 20 GeV and values shifted 
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lifetime stau [sec] gravitino mass [GcV] 



Figure 2: a) Scaled momentum spectrum p/m — of f production with contributions from fifi (magenta), eRCR^L 
(cyan), fiRp,R (yellow) and XiXi (gr^en); b) r recoil spectrum of the decay fi — > tG with Ejet of r ^ pvr + SnVr 
compared with simulations assuming itiq = 20GeV (red histogram), lOGeV (blue) and 30GeV (magenta); c) f lifetime 
distribution; d) f lifetime versus G mass for various SuSY scenarios representing different f masses. Statistics correspond 
to GDM e scenario, £ = lOOfb"^ at = 500 GeV 

by ±10 GcV arc shown as well. A fit to the 'jet' energy spectrum, using either an analytical formula for the 
endpoint or a complete simulation, yields a gravitino mass 

= 20 ± 4 GcV . (7) 

Taking all results together one can test the gravitational coupling of the stau to the gravitino and access the 
Planck scale, respectively Newton's constant. Inserting the expected values and accuracies on ?7if , tf and ttiq 
in eq. (3) one finds for the supergravity Planck scale 

Mp = (2.4 ±0.5) • 10^^ GeV , 

where the error is dominated by the precision on the gravitino mass measurement. 

The gravitino mass can be deduced more precisely from the f mass and lifetime, if the gravitational coupling 
is assumed and the macroscopic value of Mp is used in the decay- width of eq. (2) . The resulting gravitino mass 
is rriQ = 20 ± 0.2 GeV, where the error is dominated by the lifetime measurement. 
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It is a unique feature of gravitino LSP scenarios that the Planck scale can be directly measured in microscopic 
particle experiments by studying the properties of the NLSP and its decay. A further interesting test to reveal 
the nature of the gravitino as the supcrsymmctric partner of the graviton would be to determine the spin. 
This is in principle possible by studying correlations in the radiative decay f t^G [3], see also sect. 6. 
Experimentally this is quite challenging, because it requires to distinguish single photons from tt*' decays, and 
in addition the expected rates are lower by two orders of magnitude. 



5.2. mSUGRA scenario GDM C 



Due to the heavy sleptons in model GDM ^ a higher centre of mass energy of \/s = 800 GeV and an 
integrated luminosity of £ = 1 ab^^ is chosen. The /37 distribution of the produced f's (inclusive cross section 
csuSY = 5 fb) is shown in figure 3 a. It exhibits the typical behaviour of slepton production close to threshold: 
a pronounced peak at /37 ~ 0.6 from fifi pairs and another enhancement around 0.3 from br and jln cascade 
decays, resulting in 1350 (850) trapped f's in the hadron calorimeter (yoke). The analyses of the momentum 
spectrum yields a f mass of mf^ = 340.2 ± 0.2 GeV. From the decay time distribution one gets the lifetime 
tf = (1.8 ± 0.06) • 10^ s. The r jet energy spectrum of the f decay, shown in figure 3 b, extends over a wide 
range. The edge of the endpoint energy, being less pronounced than in the previous example, can be used to get 
a direct gravitino mass measurement of tuq = 100 ± 10 GeV. This value may be compared with the gravitino 
mass TOg = 100 ± 2 GeV calculated from the stau mass and lifetime. 

Note: If higher energies were accessible at the ILC, then data taking at y/s = 900 GeV, i.e. just above XiXi 
threshold, increases the number of trapped f's by a factor of 1.5 in the GDM C scenario. 




Figure 3: a) Scaled momentum spectrum p/m = Pj of f production with contributions from fifi (magenta), CReu (cyan) 
and fiRfiR (yellow) ; b) r recoil spectrum of the decay fi tG with Ej^t of r ^ pUr + SnuT compared with simulations 
using 771(5 = 100 GeV (red histogram), 80 GeV (blue) and 120 GeV (magenta). GDM C scenario, C = 1000 fb"^ at 
= 800 GeV 



The GDM rj scenario has an almost identical sparticle spectrum and differs essentially only by the lower 
gravitino mass of = 20 GeV. This causes a much shorter lifetime, but, more importantly, the r recoil energy 
spectrum is not really sensitive to the gravitino mass (mg/mf = 0.06) and can only be used to set an upper 
limit. 
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5.3. mSUGRA scenario FS 600 

The model FS 600 is investigated at a cms energy -^/s = 500 GeV, close to threshold pair production of the 
moderately heavy sleptons and neutralinos. A larger integrated luminosity is needed in order to achieve similar 
accuracies as for the other mSUGRA benchmarks. The results of the analysis are given in table III. It is 
worth noting that operating the ILC at a slightly higher energy of 520 GcV leads to higher cross sections and 
a significant increase by a factor of 1.5 for the rate of trapped t's. 

5.4. Gauge mediated symmetry brealting scenario SPS 7 

The SPS 7 scenario has a relatively light sparticle spectrum and is investigated assuming = 410 GeV and 
C = 100 fb^^, with a large cross section of asusY = 420 fb for inclusive f production. As can be seen in the 
distribution of figure 4 a, most f's leave the detector. There is, however, a large signal at ~ 0.4 from ejiCL 
production just above threshold, which contributes to the sample of 10000 (4900) trapped f's in the calorimeter 
(yoke). The analysis of the f momentum spectrum yields m^^ = 124.3 ± 0.1 GeV and from a fit to the decay 
time distribution one obtains tf ~ 209.3 ± 2.4 s. These values can be used to derive a very accurate gravitino 
mass of TTig ~ 100 ± 1 MeV assuming a gravitational coupling, see figure 2 d. The r recoil energy spectrum, 
shown in figure 4 b, is not sensitive to such low gravitino masses and can only serve to give an upper limit on 
a direct measurement of < 9 GeV at 95% confidence level. 
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Figure 4: a) Scaled momentum spectrum p/m — of f production with contributions from fifi (magenta), eRCR^L 
(cyan), fiRp,R (yellow) and XiXi (green); b) r recoil spectrum of the decay fi tG with Ej^t of r ^ pVr + ZtvVt 
compared with simulations using uiq = GeV (red histogram) and 10 GeV (blue). SPS 7 scenario, C = 100 fb"^ at 
JTs = 410 GeV 



In the present case one does not profit from the calorimetric information to asses the gravitino mass directly 
and relies entirely on the stau properties. One may afford to run at higher energy, thereby shifting the /?7 
spectrum of figure 4 a to larger values. For example, running the ILC at the 'canonical' energy y^s = 500 GeV 
reduces the fifi rate by 10% and the sample of trapped f's by a factor of two, degrading the precision on the 
f lifetime to 2 %. 

Note: The relative precision on the f lifetime measurement does not depend on the gravitino mass, should 
it be much lighter as often assumed in gauge mediated supersymmetry models. Technically, there may be a 
limitation to measure very short lifetimes below a millisecond. In the present method the f decay is required 
to occur uncorrelated to any beam collision, which essentially means outside a time interval of Ai = 1 ms for 
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a whole bunch train, repeating at a rate of 5 Hz. As an illustration, in the SPS 7 scenario a f lifetime of Sms 
corresponds to a gravitino mass of 0.5 MeV. 



5.5. Gaugino mediated symmetry breaking scenarios P3 

The model P3 is investigated at ^/s = 500 GeV and assuming £ = lOOfb"^ The inclusive f cross section of 
f^susY = 470 fb provides large rates, as can be seen in the spectrum of figure 5 a. The main contributions 
to the sample of stopped f's come from XiXi f-^d srcl production, one expects about 3900 and 3700 particles 
in the calorimeter and yoke, respectively. From the kinematics of fifi production one gets a f mass of m^^ = 
102.5 ± 0.2 GeV. A fit to the f lifetime distribution yields tf = (4.2 ± 0.1) • 10^ s. A direct measurement of 
the gravitino mass using the r recoil energy spectrum of f ^ tG, as displayed in figure 5 b, provides a value 
of = 25 ± 1.5 GeV. Alternatively, one derives a gravitino mass of m^, = 25 ± 0.3 GeV from the f lifetime 
measurement. 

Note: The number of trapped f's in the calorimeter can be increased by a factor of 1.8 when running at a 
slightly lower cms energy y/s = 480 GeV. 




Figure 5: a) Scaled momentum spectrum p/m — /37 of f production with contributions from fifi (magenta), enfiR^L 
(cyan), jlRfiR (yellow) and XiXi (green); b) r recoil spectra of the decay fi tG with Ejet of t — > + 3tvUt compared 
with simulations assuming mg, — 25 GeV (red histogram), 20 GeV (blue) and 30 GeV (magenta). x^SB P3 scenario, 
£ = 100 fb"^ at ^ = 500 GeV 

The scenario P3 yields a lifetime of about 1.3 years. Taking instead a larger gravitino mass of m^, ~ 50 GeV, 
as proposed by the authors of [14], the lifetime increases by almost an order of magnitude, and it may be 
difficult to observe all accumulated f's during the operation of an experiment. Assuming a measurement period 
of three years will provide a lifetime of tf ~ (3.9 ± 0.7) • 10® s, corresponding to 12.4 years. The relative error 
grows faster than expected from the reduced statistics (factor 4 less decays) when observing only a fraction of 
the lifetime. Of coarse, this has also consequences for the measurement of the decay spectra. 

The xMSB PI scenario is characterised by somewhat heavier sparticles and is investigated at y/s = 420 GeV, 
just above threshold production of sleptons and the lightest neutralino. Again, the event statistics allow mea- 
surements to be performed at the per mill to per cent level, see table III. 



11 



6. Discussion of results 

If kinematically accessible, metastable staus are copiously produced at the ILC. They are easy to detect 
and their properties, like mass and lifetime, and the decays into gravitinos can be studied in great detail. The 
results of the various case studies of the supersymmetry models and the assumed experimental conditions are 
summarised in table III. The common features of all spectra concerning the achievable accuracies of the basic 
measurements are: (i) the f mass can be measured at the per mill level by exploiting the 2-body kinematics 
of fifi pair production in the tracking detector; (ii) the f lifetime can be determined within a few per cent 
by observing the decays f tG triggered in the calorimeter and/or iron yoke; and (iii) the G mass can be 
accessed directly with a precision of order 10 % (degrading rapidly for m^/mf < 0.1) by measuring the r recoil 
spectra of the f decays in the calorimeter. 

Table III: Expected accuracies on the determination of f and G properties for various SuSY scenarios: stau mass mf, 
stau lifetime tf, gravitino mass ni?^ from r recoil energy spectra and gravitino mass deduced from f lifetime. 
Experimental conditions in the last columns: centre of mass energy ^/s, inclusive f production cross section (if fx, 
integrated luminosity C, number of trapped t s in hadron calorimeter 7V|"^^' and instrumented yoke NY"^" 





mf [GeV] 


tf [s] 


m^- [GeV] 


m'T [GeV] 


[GeV] 


o-f f X [fb] 


£[fb-i] 




f 


1 SPS 7 


124.3 ±0.1 


209.3 ± 2.4 


< 9 


0.1 ±0.001 


410 


420 


100 


10000 


4900 


2 FS 600 


219.3 ±0.2 


(3.6 ±0.1) 10'' 


50 ±9 


50 ±0.7 


500 


20 


250 


2100 


4200 


3 GDM e 


157.6 ±0.2 


(2.6 ±0.05) lO'^ 


20 ±4 


20 ±0.2 


500 


300 


100 


4100 


1850 


4 GDM C 


340.2 ± 0.2 


(1.8 ±0.06) 10® 


100 ± 10 


100 ±2 


800 


5 


1000 


1350 


800 


5 GDM 77 


322.1 ±0.2 


(6.9 ±0.3) 10'' 


20 ±25 


20 ± 0.4 


800 


9 


1000 


1050 


1850 


6 xMSB PI 


185.2 ±0.1 


(9.1 ±0.2) 10'' 


50 ±3 


50 ±0.6 


420 


27 


200 


3700 


4100 


7 xMSB P3 


102.5 ±0.2 


(4.2 ±0.1) lO'^ 


25 ± 1.5 


25 ±0.3 


500 


470 


100 


3900 


3700 



These expectations arc based on moderate integrated luminosities, the event rates to be accumulated during 
one to three years of data taking at the ILC under nominal conditions. Although the f decays are treated in 
a parametrised form the results should be reliable within a factor of two. A detailed detector simulation may 
lead to a reduction of the data samples due to acceptance cuts, quality criteria and resolution effects. However, 
the losses may be compensated and the event statistics can be considerably increased by optimising the beam 
energies (see sect. 5) and making use of the beam polarisations. For instance, the production cross sections for 
light sleptons, fifi, eReR and /ifl/ifl, are enhanced by a factor (l + 'Pg-)(l — 7'g+) ~ 2.9 when using right-handed 
polarised electrons of degree V^- = +0.8 and left-handed polarised positrons of Ve+ = —0.6. 

Concerning the expected ILC detector performance with respect to detecting metastable f's there is, however, 
a caveat. At present, it is envisaged to operate the LDC calorimeters in a pulsed mode [19]. In order to cope 
with the heat production of the electronics, it is foreseen that the amplifiers and read out channels of the 
detector parts will be switched on and read out just during one bunch train for 2 ms, and then be switched off 
to wait for the next bunch train occurring 200 ms later; i.e. the detector is inactive for most of the time. Such 
an operation is designed for usual, beam correlated physics scenarios, but clearly has to be revised if one wants 
to observe the decays of metastable particles. For the iron yoke instrumentation it does not really seem to be 
necessary to apply a power pulsing. A quasi-continuous read-out of all channels above a certain threshold over 
the so far 'idle' time gap of 198 ms appears possible, replacing the f decay trigger by software analysis [26]. 
Thus a lifetime measurement will be feasible in any case, may be at the expense of a lower event rate restricted 
to f's trapped in the yoke. The goal to also increase the duty cycle of the hadron calorimeter substantially 
has to be pursued in future R&D detector developments. Since a direct measurement of the gravitino mass is 
extremely important, it will be assumed in the following discussion that the technical problems will be solved 
and the calorimeter will be permanently sensitive with essentially no dead-time. 
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An independent measurement of the gravitino mass is of prime importance. Applying eq. (3) one can deter- 
mine the supergravity Planck scale from the NLSP decay with an accuracy of order 10%. A value consistent 
with the macroscopic scale of gravity, i.e. Newton's constant Gjv = 6.7 • 10~^^GcV~^, would be a decisive 
test that the gravitino is indeed the supcrpartncr of the graviton. At the same time the measurement of the 
gravitino mass can be used to yield the supersymmetry breaking scale F = \/3A/p 7713/2, see eq. (1), which is 
an important parameter to unravel the nature of the supersymmetry breaking mechanism. 

From the f mass and lifetime measurements one can compute the gravitino mass assuming that the f decay 
occurs with gravitational coupling according to eq. (2). The expected accuracy is again of order per cent. 
This interpretation may be plausible, however, it can only be justified, if there is additional information which 
confirms the gravitino LSP assignment. 

For too light gravitinos a direct mass measurement may not be possible, leaving the nature of the LSP open. 
An alternative option may be a spin 1/2 axino d as lightest supersymmetric particle [27]. In such a scenario the 
f lifetime is not related to the axino mass and, depending on the parameters, may vary between some 0.01 sec 
and 10 hours, comparable to the expectations of gauge mediated and gravity mediated SuSY breaking models. 
A method to differentiate between these two scenarios is to study the radiative three-body decays f — > tjG 
and f — > Tja. The branching ratios are suppressed by two orders of magnitude and expected to be slightly 
larger for the axino LSP, the exact values depend on the selected phase space to identify an isolated photon. 
More promising is to investigate 7 — r correlations of the decay topology. In the case of gravitino LSP the 
photons are preferentially emitted coUinear with the r and have a soft, Bremsstrahlung like energy spectrum. A 
similar configuration also occurs for the axino LSP, however, there is in addition a substantial rate of energetic 
photons, E-y —>■ mf/2, emitted opposite, back-to-back to the r direction. Such a signature would clearly allow 
to distinguish between the interpretations as axino LSP or as gravitino LSP. Experimentally the analysis of 
radiative f decays is quite ambitious because one has to discriminate a single photon against the photons from 
TT*^ decays as well as against the hadrons in hadronic r decays. The expected low event rates will be further 
reduced by efficient selection criteria. Hence, large statistics data samples are required, which may be provided 
by the ILC. 

7. Conclusions 

It has been shown in detailed analyses that future ILC experiments have a rich potential to study SuSY 
scenarios where the gravitino G is the lightest supersymmetric particle and a charged stau f is the long-lived, 
metastable next-to-lightest supersymmetric particle. Precise determinations of the f mass and lifetime and 
of the G mass appear feasible, provided the proposed detectors will operate at a reasonable duty cycle. A 
calorimetric measurement of the gravitino mass from the r recoil spectra of the decay f —^ tG gives access to 
the gravitational coupling, i.e. to the Planck scale, in a microscopic particle experiment and thus provides a 
unique test of supergravity. Furthermore these observations will put stringent constraints on an interpretation 
of the gravitino as dark matter candidate, being undetectable in astrophysical search experiments, and will 
allow its relic density to be to computed reliably. 

Acknowledgements. I want to thank W. Buchmiiller and P. M. Zerwas for many valuable discussions and 
comments on the manuscript. 

References 



[1] H. Pagels and J. R. Primack, Phys. Rev. Lett. 48 (1982) 223; M. Y. Khlopov and A. D. Linde, Phys. Lett. 
B 138 (1984) 265; J. R. Ellis, J. E. Kim and D. V. Nanopoulos, Phys. Lett. B 145 (1984) 181; M. Bolz, 
W. Buchmiiller and M. Pliimacher, Phys. Lett. B 443 (1998) 209 [arXiv:hep-ph/9809381]. 



13 



[2] S. Eidelman et al. [Particle Data Group], Phys. Lett. B 592 (2004) 1. 

[3] W. Buchmiiller, K. Hamaguchi, M. Ratz and T. Yanagida, Phys. Lett. B 588 (2004) 90 [arXivdiep- 
ph/0402179]. 

[4] M. Bolz. A. Brandenburg and W. BuchmuUcr, "Thermal produetion of gravitinos," Nucl. Phys. B 606 

(2001) 518 [arXiv:hep-ph/0012052]. 
[5] J. L. Feng, A. Rajaraman and F. Takayama, Phys. Rev. Lett. 91 (2003) 011302 [arXiv:hep-ph/0302215] 

and Phys. Rev. D 68 (2003) 063504 [arXiv:hep-ph/0306024]. 
[6] K. Hamaguchi, Y. Kuno, T. Nakaya and M. M. Nojiri, Phys. Rev. D 70 (2004) 115007 [arXivdiep- 

ph/0409248]. 

[7] J. L. Feng and B. T. Smith, Phys. Rev. D 71 (2005) 015004 [Erratum-ibid. D 71 (2005) 0109904] [arXivdiep- 
ph/0409278]. 

[8] A. De Roeck, J. R. Ellis, F. Gianotti, F. Moortgat, K. A. Olive and L. Pape, arXiv:hep-ph/0508198. 

[9] M. Dine, A. E. Nelson and Y. Shirman, Phys. Rev. D 51 (1995) 1362 [arXiv:hep-ph/9408384]. 
[10] T. Moroi, H. Murayama and M. Yamaguchi, Phys. Lett. B 303 (1993) 289. 
[11] B.C. AUanach et al, Eur. Phys. J. C 25 (2002) 113 [arXiv:hep-ph/0202233]. 
[12] H. P. Nilles, Phys. Rept. 110 (1984) 1. 

[13] D. E. Kaplan, G. D. Kribs and M. Schmahz, Phys. Rev. D 62 (2000) 035010 [arXiv:hep-ph/9911293]; 

Z. Chacko, M. A. Luty, A. E. Nelson and E. Ponton, JHEP 0001 (2000) 003 [arXiv:hep-ph/9911323]. 
[14] W. Buchmiiller, J. Kersten and K. Schmidt-Hoberg, JHEP 0602 (2006) 069 [arXiv:hep-ph/0512152]. 
[15] W. Buchmiiller, K. Hamaguchi and J. Kersten, Phys. Lett. B 632 (2006) 366 [arXiv:hep-ph/0506105]. 
[16] W. Porod, Comput. Phys. Commun. 153 (2003) 275-315 [arXiv:hcp-ph/0301101]. 
[17] A. Djouadi, J. L. Kneur and G. Mouhaka, arXiv:hep-ph/0211331. 
[18] Tesla Technical Design Report, DESY 2001-011, Part IV: A Detector for TESLA. 
[19] Large Detector Concept working group, http://www.ilcdc.org 
[20] B Rossi, 'High-Energy Particles', Prentice-Hall, Inc (1952). 

[21] T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, S. Mrenna and E. Norrbin, Comput. Phys. 

Commun. 135 (2001) 238 [arXiv:hep-ph/0010017]. 
[22] T. Ohl, Comput. Phys. Commun. 101 (1997) 269 [arXiv:hep-ph/9607454]. 
[23] M. Pohl and H. J. Schreiber, DESY-02-061, arXiv:hep-ex/0206009. 

[24] A. Freitas, A. von Manteuffel and P. M. Zerwas, Eur. Phys. J. C 34 (2004) 487 [arXiv:hep-ph/0310182], 
and Eur. Phys. J. C 40 (2005) 435 [arXiv:hep-ph/0408341]; A. Freitas, D. J. Miller and P. M. Zerwas, Eur. 
Phys. J. C 21 (2001) 361 [arXiv:hep-ph/0106198]. 

[25] J. L. Feng and M. E. Peskin, Phys. Rev. D 64 (2001) 115002 [arXiv:hep-ph/0105100]. 

[26] G. Eckerlin, private communication. 

[27] A. Brandenburg, L. Covi, K. Hamaguchi, L. Roszkowski and F. D. StefFen, Phys. Lett. B 617 (2005) 99 
[arXiv:hep-ph/0501287] . 



